Propagation and Relaxation of Tension in Stiff Polymers 
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We present a unified theory for tfie longitudinal dynamic response of a stiff polymer in solution to 
various external perturbations (mechanical excitations, hydrodynamic flows, electrical fields, tem- 
perature quenches . . . ) that can be represented as sudden changes of ambient /boundary conditions. 
The theory relies on a comprehensive analysis of the non-equilibrium propagation and relaxation of 
backbone stresses in a wormlike chain. We recover and substantially extend previous results based 
on heuristic arguments. New experimental implications are pointed out. (LMU-ASC 15/05) 
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Despite considerable practical and interdisciplinary in- 
terest, it is theoretically not yet fully understood how 
polymers respond to external fields Consider, e.g., 

the simple problem of an inextensible semiflexible poly- 
mer suddenly stretched along its end-to-end vector by an 
external force f [Pulling). It was recognized before Q 
that tension propagation (from the ends into the bulk) is 
the key to understanding its dynamics: in response to the 
spreading tension, the polymer stretches within a grow- 
ing boundary layer of length £|| , as depicted in Fig. ^ 
Depending on the setup, different tension propagation 
laws £||(i) have been predicted 0, 0, S 0- particu- 
lar, we contrast the above Pulling-Bccwaiio with the (in- 
verse) i?e/eas e-scenario, where a constant pre-strctching 
force f is suddenly removed. While for small f, one ex- 
pects i\\{t) (X t^^^ in both cases 0, the predictions for 
strong force are markedly different: £\\{t) oc {ft)^^^ 
for Pulling, and ^ (t) cx f /^tVa [| 

for Release. However, 
the precise meaning of "strong" and "weak" and the va- 
lidity of the diverse assumptions in Refs. 0, 0] are not 
obvious. Here we develop from first principles a theory of 
stress propagation and relaxation that allows us to unify 
and systematically extend these previous results, and to 
derive (often analytically) the longitudinal nonlinear re- 
sponse to various external perturbations. After intro- 
ducing the standard model of a semiflexible polymer, we 
extend a heuristic argument of Ref. Q for Pulling. This 
elucidates the crossover from "weak-" to "strong-force" 
behavior and reveals the crucial length-scale separation 
underlying our subsequent systematic analysis. 

The wormlike chain model represents the polymer by a 
smooth inextensible contour r(s, t) subject to the energy 
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The scalar force /(s,t) (backbone "tension") is a La- 
grange multiplier function introduced to enforce the 
inextensibility constraint r'^ = 1 for the tangent vec- 
tor r' = dr/ds. We require a bending stiffness k such 
that the persistence length £p = K/ksT is much larger 




FIG. 1: Pulling (schematic): In response to an external force 
f, the thermally undulated contour r(s) — {rj_,s — ry)"^ is 
straightened within boundary layers of growing width £\\{t). 



than the contour length L, which entails relative mean 
square transverse displacements of order e = L/£p <^ 1. 
The elastic forces derived from H have to balance ther- 
mal forces ^ (Gaussian white noise) and Stokes fric- 
tion, which (in the "free-draining" approximation) enters 
through two local friction coefficients per unit length , 
Cjt ~ for motion perpendicular and parallel to r', 

respectively: [Cy'''''''' + C±(l ~ r'r')] ■ dtr = —STC/Sr + ^. 
For the following, we choose convenient units such that 
K = = 1. Then, all dimensional quantities represent 
powers of length (e.g. ksT = £~^) and = sa 1/2. 

We now turn to a heuristic discussion of Pulling to 
leading order in e. In the parameterization introduced 
in Fig. ^ the exact equations of motion reduce to an 
equation for the transverse displacements r± alone, 
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with a spatially uniform tension / ~ f B(<) fixed by the 
driving force at the boundaries; {i±{s,t) ■ ^±{s',t')) = 
4£~^S{s ~ s')S{t — t'). The (higher order) longitudinal 
displacements r|| are slaved by the arclength constraint 
r||=rf/2. 

From a simple scaling analysis of Eq. (0), {'f±)/t w 
{r±_) {£^'^ + f ^_i^)7 we deduce the characteristic dynamic 
wavelength £±{t) corresponding to the (lowest) mode 
equilibrated at time t. For instance, £±{tj;) = L defines 
the longest relaxation time. Due to the competition be- 
tween bending forces (oc r±£j^'^) and tension (cx rj_f ^^^), 
the growth of £± exhibits a dynamic crossover from 
free relaxation {£± oc i^/^) to relaxation under tension 
{£± oc \/ft) at a characteristic time if = (Tab.p/left). 
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TABLE I: The transverse equilibration length £±{t) and the 
tension propagation length £\\ (t) both exhibit a crossover at 
tf = (here, for the Pulling problem with f ^ L~^, t <^tj;). 
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By the above interpretation of £±, the longitudinal 
elongation of a subsection of arclength £± under a given 
tension f can be estimated by equilibrium theory. One 
has to distinguish weak and strong tension relative to 
the internal characteristic force scale ^J^^ of the subsec- 
tion, which corresponds to the Euler buckling force of 
the subsection. For weak tension f ^J7^! the elongation 
f£'^£~^ follows from linear response |8j. For strong ten- 
sion f ^ f|^^ the subsection is virtually straight, so that 
the elongation is equal to its equilibrium thermal contrac- 
tion £'j_£p^ caused by the bending undulations. Since the 
whole polymer is subject to the same constant tension f, 
it can be divided (at any time t) into L/£±(t) indepen- 
dent equilibrated subsections of length £±{t). The total 
change AR{t) = \{R{t) - R{0))\ of the end-to-end dis- 
tance R thus crosses over from AR oc Lf£'[£~^ oc t^^^ 
for t < if to Ai? oc L£^£p'^ oc t^/^ for t > tf. 

These results comprise the predictions of ordinary per- 
turbation theory (OPT) to leading order. As evident 
from Eq. (0), it neglects longitudinal friction forces, 
which arc of higher order in e. However, the re- 
sulting AR obtained above implies a total longitudi- 
nal friction C,LdtAR crossing over from Ci^fi'^/^C^ to 
CL2(f/t)i/2^-i at tf. Both expressions diverge for 
t ^ Q indicating the breakdown of OPT at short times. 
More precisely, for times shorter than a certain ti, 
the longitudinal friction would exceed the driving force f. 
This was recognized as a consequence of the omission of 
tension propagation in Eq. |5J): it was argued IJI that ac- 
tually only segments up to a distance £\\ (t) from the 
ends are set into longitudinal motion causing longitudinal 
friction. The proper expression for the total longitudinal 
friction thus follows from the above upon substituting L 
by £|| . The balance of longitudinal friction and external 
force can now be satisfied by choosing the size £\\{t) of 
the boundary layer according to Tab.p/right. Hence, the 
putative "weak- and strong- force" cases £« oc i^/^ j3| and 
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^11 oc t^/^ jsj are identified as asymptotes of a "short-long 
time" crossover (still) signaling the change from "free" 
to "forced" relaxation at t ^ tf. 

In summary, the foregoing discussion reveals: (i) the 
breakdown of OPT at times t < t^, where (ii) the ne- 
glected longitudinal friction actually limits the relaxation 
to boundary layers of size £|| ; (Hi) the crossover from free 
to forced relaxation at t = tf; (iv) the scale separation 
i±/£\\ oc < 1 (Tab.©. 



The latter lends itself as starting point for a multiple- 
scale calculus to separate the physics on different dy- 
namic scales and obtain an improved ("multiple-scale") 
perturbation theory (MSPT) that is regular in the limit 
i — > while e <C 1 is fixed. The procedure is similar to 
that for athermal dynamics and will be documented 
in detail elsewhere [l^l . The basic idea is to regard func- 
tions g{s) as functions g{s,se^^^) of rapidly and slowly 
varying arclength parameters s and se-^^^ that are treated 
as independent variables. Closed equations for the dy- 
namics on the scale se^/^ are obtained upon averaging 
g{s) = J^ds g{s,se^^^)/l over the fluctuations on the mi- 
croscale (£± 'ti I 't^ £\\). To leading order, we get 

f = mt), 



dtr^ = -r'l" + fr'l + , and dif = -(dtr^ . (3) 

This provides the sought-after rigorous local general- 
ization of the above, heuristically renormalized force 
balance. The arclength average extending over many 
{l/£± ^ 1) uncorrelated sections of length £± subject 
to the same tension /(s) can be interpreted as a coarse- 
graining that effectively generates a local ensemble av- 
erage: 9tr|| ^ (dtf'^^) for e 0. Only the "system- 
atic" 0(1)— variations of the tension are retained, while 
its ©(e)— noise is leveled out, so that the longitudinal 
Eq. (jSl represents deterministic dynamics: local longi- 
tudinal motion is driven by tension gradients (like in a 
thread pulled through a viscous medium) . 

Integrating the longitudinal Eq. Q over time expresses 
the change of the thermal contraction ry in terms of the 
time-integrated tension F = J^dt' f, namely (^r'^^{s,t) — 
rj|(s,0)) = -dlF/C. Since 2r|j = rf from the arc- 
length constraint, we integrate the transverse Eq. ^ for 
{r'l){-s,t) = {[^Tj'^Jt' qX^{T,t,t')iAl,t')f)- Here 



xx{q;t,t') 



(4) 



should be recognized as the susceptibility for the response 
of the Fourier modes of r±_ to transverse forces. Note 
that the s— dependence of (r^) is purely adiabatic, as 
it is parametrically inherited from F(s,t). Alto get her, 
Eq. (j2J) is condensed into a single equation for F [l8j , 
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dt' xl{q;t,t') 



(5) 

Indices "— " referring to t < were introduced to al- 
low the system to be prepared in equilibrium with ambi- 
ent/boundary conditions different from those prescribed 
for t > 0. Taking £p (for t > 0) to oo, only the first term 
in the integrand remains, which thus accounts for the de- 
terministic relaxation of the initial thermal contraction 
(set by £~ , f~). In this limit, the "zero-temperature" 
buckling dynamics analyzed in Ref. is recovered. For 
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finite £p the second term in the brackets represents the 
action of the thermal forces for t > 0, which aim to es- 
tablish the equilibrium contraction. 

To further unravel the physical implications of Eq. ij^l, 
we begin with the tension propagation regime ^|| <^ L, 
where L is irrelevant. Problems like Pulling and Release 
still depend on four length scales {i~ = £p, f~^/^, s, t^/^). 
Yet, Eq. jSJ is solved exactly by the two-variable scaling 
ansatz 

/(s,t)=f<i>(sM,t/tf) . (6) 

With = f-2 and Sf = {lp/C,f''^t]'^ Eq. ® can be 
shown to eliminate the parameter dependencies in Eq. © 
and in its boundary conditions. The scaling function $ 
describes how sudden changes of the tension at the ends 
spread into the bulk of the polymer. In the limits t <^t^ 
and < » Eq. © reduces to the simple (one-variable) 
scaling laws $ ~ {t/t^)°'(j){s/t\\) with the tension prop- 
agation length = sj(t/tj)^. Note that this crossover 
scenario, as well as the expressions for tj, sj w ^||(if), and 
^11 are consistent with our above heuristic observations for 
Pulling (Tab.Hl). In fact, this structure generally emerges 
for all problems involving a single external force scale. It 
is implicitly understood that (/>, a, and z will generally 
not only depend on the kind of external perturbation 
applied, but will also be different in both scaling limits 
i/if 0,00. In the following, these limits are analyzed 
in more detail. 

For t Eq. (jSJ may be linearized in / and the 

scaling function can be obtained analytically [l^ . In 
contrast to (j), the corresponding exponent z = 1/8 does 
not depend on the boundary conditions. It already fol- 
lows from requiring (j) to become f— independent, as in 
linear response. The short-time dynamics for strong ex- 
ternal force is thus closely related to the linear response. 
As established by our heuristic discussion of Pulling^ this 
is due to the relaxation of subsections with Euler forces 
fj^^ much larger than the external force. Note, however, 
that the limit f is problematic, as it does not in- 
terchange with e — > [ijj- In fact, extending Eq. 10) to 
linear response amounts to an uncontrolled factorization 
approximation {fr]_) (/) (r^)- 

For t ^ tf the dynamics becomes nonlinear in the 
external force and starts to depend on the force proto- 
col. Previously predicted power laws can be recovered 
from Eq. ij^l by employing different approximations to 
its right hand side. In the taut-string approximation of 
Ref. |3| one neglects for t > bending and thermal forces 
against the tension, i.e. one drops the g^— term in the 
expression Eq. (0J for x± and sets ip ^ oo. The comple- 
mentary quasi-static approximation of Ref. 0| amounts 
to the omission of memory effects, i.e. to the assump- 
tion of instantaneous equilibration of tension and stored 
length, F{-t}_- F(t') f(t){t - t'). A careful analysis 
of Eq. jSJ |l3l shows that either of these approximations 
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FIG. 2: Double-logarithmic sketch of the tension propagation 
laws l\\{t) (X . At tf = f^^ they cross over from a univer- 
sal short-time regime to (problem-specific) tension-dominated 
intermediate asymptotics, except for weak forces, f < l^/L'^, 
and for Ip — Quenches. The propagation ends when l\\{t) ^ L. 



becomes rigorous in the intermediate asym pto tic regime 
defined by t > tf, £\\ < L. The "pure" ^ scenarios 
are summarized in Fig. [21 and for the cases Pulling and 
Release also in Fig. 13 and Tab. ITTl 

In Fig. 13 we have moreover displayed results for sud- 
den changes in persistence length from £~ to £p ^ £~ 
{"£p— Quench") and Electrophoretic Pulling, which have 
not been discussed before. The second is a variant of 
the Pulling problem, where the external force is applied 
along the whole contour of an end-grafted polymer as is 
the case for a hydrodynamic flow or an electric field. The 
£p- Quench is exceptional in that there is no external force 
scale involved, so that Eq. (|SJ) can be solved by a simple 
one-variable scaling ansatz f{s,t) = t~^/^ip{s/£\^) with 
£|| « {£p/(^y/^t^^^ . Neither the taut-string approxima- 
tion nor the quasi-static approximation can be applied, 
and in contrast to the other cases the scaling function 
has to be evaluated numerically. 

Eventually, at a time the tension will have propa- 
gated through the whole polymer, i.e. ^[[(il) ~ L. Subse- 
quently, the characteristic longitudinal scale is the con- 
tour length L. One would expect that regular pertur- 
bation theory would then become valid. Surprisingly, 
for ^ if the Release scenario provides an excep- 
tion. The contraction dynamics exhibits an intermedi- 
ate regime of homogeneous tension relaxation (white in 
Fig. |3J). Its asymptotic power-law form is revealed by 
the same quasi-static approximation applicable during 
the preceding tension propagation, but with the sepa- 
ration ansatz f{s,t) ~ h{s){(L'^/£pt)'^^^ instead of the 
single-variable scaling ansatz. It solves Eq. (jsj analyti- 
cally with a roughly parabolic stationary tension profile 
h{s). The homogeneous tension relaxation dominates the 
short-time relaxation for t <^ ti, = {(^L?' /£pY if f ^ cxd 
(i.e. <i ^ 0). 

To make contact with experiments, it is desirable to 
derive the consequences for pertinent observables such 
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Release ^ 5 Pullinq 

FIG. 3: Characteristic times (logarithmic scale) for Pulling 
and Release against the applied external force f (linear scale). 
The time (stars) separates regions where ordinary pertur- 
bation theory (OPT) applies (dark shaded) from regions (light 
shaded) of linear (hatched) and nonlinear tension propagation 
and from homogeneous tension relaxation (white). Whereas 
longitudinal friction is negligible for t > t*, it limits the dy- 
namics for t <ti,. 



TABLE II: Growth laws for the end-to-end distance A_R(t) 
in the intermediate asymptotic regimes marked in Fig. |3] 
OPT and MSPT refer to "ordinary perturbation theory" and 
"multiple-scale perturbation theory" , respectively. 





Release 1 Pulling 


linear MSPT 
linear OPT [g] 
nonlin. MSPT 
hom. relaxation 
nonlin. OPT 


2ir(f)-i( 
2-fr(2r 

3.503 (C4)"M^i^ 
2-54S(L/C^p)3t3 


C^p)-Mt5 

"i(L/^p)ftf 

(2/^)5 (L/^,)(ff)^ 



toskeleton |y, [iSj. Inclusion of liydrodynamic interac- 
tions merely produce logarithmic corrections but would 
give rise to more interesting effects for membranes. Other 
natural generalizations including the transverse nonlinear 
response and more complex force protocols (e.g. [20]) are 
currently also under investigation. 



as the end-to-end distance. Integrating the longitudi- 
nal Eq. (0) over s and t shows that the growth laws 
A^(t) = 2C~^ \F'{s = 0,t)| are directly related to the 
tension profiles discussed above. Tab. IITI summarizes our 
results for the intermediate asymptotic regimes. Note 
that is a coarse-grained quantity that does not re- 
solve the "microscopic" details below the coarse-graining 
scale I. Near the polymer ends these give relevant con- 
tributions obliterating the predicted t^^^ in experiments 
0, 01 . During homogeneous tension relaxation A.R oc 
t^/^, which we expect to hold for strongly stretched poly- 
mers even if L ^ f p (e.g. DNA), at variance with earlier 
predictions The exponent 1/3 coincides with that 
obtained by adiabatic application of the stationary force- 
extension relation to a "frictionless" polymer with 
attached beads at its ends. Finally, AE « (C^p)"^/^i^/® 
in Ip— Quenches for t ^ t'[. Interestingly, the tension 
propagation/relaxation itself can in some situations be 
directly monitored experimentally. In ip- Quenches, the 
(longitudinal) radius of gyration mirrors the character- 
istic bulk relaxation / cx t^^/^ of the tension In 
Electrophoretic- Pulling, where AR oc t for t ^ the 
force on the grafted end obeys f oc £|| . 

In conclusion, we have developed a unified theory of 
non-equilibrium tension dynamics in stiff polymers based 
on the scale separation between the two dynamic corre- 
lation lengths £± and ^||. The recovered known results 
and our new predictions are summarized in Figs. |21 El 
and Tab. |n| Various dynamic regimes should be well 
realizable for certain biopolymers and it is an intrigu- 
ing question, whether the tension propagation laws €|| (t) 
govern mechanical signal transduction through the cy- 



[1] T. T. Perkins, D. E. Smith, and S. Chu, Science 276, 
2016 (1997); B. Maier, U. Seifert, and J. O. Radler, Euro- 
phys. Lett. 60, 622 (2002); N.-K. Lee and D. Thirumalai, 
Biophys. J. 86, 2641 (2004). 

[2] Y. Bohbot-Raviv, W. Z. Zhao, M. Feingold, C. H. Wig- 
gins, and R. Granek, Phys. Rev. Lett. 92, 098101 (2004). 

[3] U. Seifert, W. Wintz, and P. Nelson, Phys. Rev. Lett. 
77, 5389 (1996). 

[4] R. Everaers et al., Phys. Rev. Lett. 82, 3717 (1999). 

[5] F. Brochard-Wyart, A. Buguin, and P. G. de Gennes, 
Europhys. Lett. 47, 171 (1999). 

[6] V. Shankar, M. Pasquah, and D. C. Morse, J. Rheol. 46, 
1111 (2002). 

[7] R. E. Goldstein and S. A. Langer, Phys. Rev. Lett. 75, 
1094 (1995). 

[8] F. MacKintosh, J. Kas, and P. Janmey, Phys. Rev. Lett. 

75, 4425 (1995). 
[9] R. Granek, J. Phys. II France 7, 1761 (1997); F. Gittes 
and F. C. MacKintosh, Phys. Rev. E 58, R1241 (1998). 

[10] D. Morse, Macromol. 31, 7044 (1998). 

[11] O. Hallatschek, E. Frey, and K. Kroy, Phys. Rev. E 70, 
031802 (2004). 

[12] O. Hallatschek et al., in preparation. 

[13] L. Le Golf et al., Phys. Rev. Lett. 89, 258101 (2002). 

[14] C. Bustamante et al.. Science 265, 1599 (1994). 

[15] M. L. Gardel et al., Phys. Rev. Lett. 91, 158302 (2003). 

[16] = {C,L^/£pf if f < e-l/L'' and U = {CL^/lpfr^ if 
f >4/L'*. See Fig.H 

[17] Identifying ^x,|| as correlation length of 9trx,||, respec- 
tively, justifies the notation. 

[18] The right-hand side of this equation has recently inde- 
pendently been derived in a related context 0. 

[19] Our identification of arclength averages with (local) en- 
semble averages after Eq. ^ breaks down for f < 
{(,/ tpf'^'^t''^ , where fiuctuations in the tension become 
comparable to its average 

[20] We postpone "mixed" scenarios (e.g. strong stretching 



5 



of a pre-stretched polymer) that involve more than one 0, , they lead to multiple crossovers even for t ^ tf. 

external force scale. Contrary to previous expectations 



